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ABSTRACT 

In papers I and II in this series, we presented hydro-dynamical simulations of jet models with 
parameters representative of the symbiotic system MWC 560. These were simulations of a pulsed, 
initially underdense jet in a high density ambient medium. Since the pulsed emission of the jet creates 
internal shocks and since the jet velocity is very high, the jet bow shock and the internal shocks are 
heated to high temperatures and should therefore emit X-ray radiation. In this paper, we investigate 
in detail the X-ray properties of the jets in our models. We have focused our study on the total 
X-ray luminosity and its temporal variability, the resulting spectra and the spatial distribution of 
the emission. Temperature and density maps from our hydrodynamical simulations with radiative 
cooling presented in the second paper are used together with emissivities calculated with the atomic 
database ATOMDB. The jets in our models show extended and variable X-ray emission which can 
be characterized as a sum of hot and warm components with temperatures that are consistent with 
observations of CH Cyg and R Aqr. The X-ray spectra of our model jets show emission line features 
which correspond to observed features in the spectra of CH Cyg. The innermost parts of our pulsed 
jets show iron line emission in the 6.4 - 6.7 keV range which may explain such emission from the central 
source in R Aqr. We conclude that MWC 560 should be detectable with Chandra or XMM-Newton, 
and such X-ray observations will provide crucial for understanding jets in symbiotic stars. 
Subject headings: circumstcllar matter — hydrodynamics — ISM: jets and outflows - binaries: sym- 
biotic - methods: numerical — X-rays: ISM 
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1. INTRODUCTION 

Highly collimated fast outflows or jets are common 
in many astrophysical objects of different sizes and 
masses: active galactic nuclei (AGN), X-ray binaries 
(XRBs), young stellar objects (YSO), pre-planetary neb- 
ulae (PPN), supersoft X-ray sources and symbiotic stars. 
In the last two objects, the jet engine consists of an ac- 
creting white dwarf. In symbiotic stars, the companion 
is a red giant undergoing strong mass loss. More than 
one hundred symbiotic stars are known, but only about 
ten systems show the presence of jets. The most famous 
systems are R Aquarii, CH Cygni, and MWC 560. 

R Aquarii, with a distance of about 200 pc, is one of the 
nearest symbiotic stars and a well known jet source. The 
jet has been extensiv ely observed in the optical and at ra- 
dio wavelengths (e .g. Solf & Ulrich 1985; Paresce fe Hack! 
11994 iHollis et all ll985al|bT ). It shows a rich morphol- 
ogy, e.g. a series of parallel features in the jet and the 
counter-jet, extending to a few hundred AU each. This 
is a hint of pulsed ejection of both jets. Furthermore, R 
Aqr is the first je t in a symbiotic system, which was de- 
tected in X-ravs (IViotti et al.lll987| iHunsch et aT1ll998l ; 
IKellogg et afl 120011 ). iKellogg et ali ((20011) found peaks 
of O VII at 0.57 keV in both the NE and SW jets and 
a peak of N VI at 0.43 keV only in the NE jet. The 
spectra are consistent with a soft component with kT ~ 
0.25 keV. The central source shows a supersoft blackbody 
emission with kT ~ 0.18 keV and a Fe Ka line at 6.4 
keV which suggests th e presence of a hard source near 
the hot star. Recently, IKellogg et al.l (|2006[ ) reported on 
five years of observations with Chandra and were able 
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to measure the pro per motion of knots in the NE jet of 
about 600 km s" 1 . iNichols et ail (|2007l ) investigated the 
X-ray emission from the inner 500 AU of this system. 

In 1984/85, CH Cygni showed a strong radio outburst, 
during which a do uble-sided jet with multiple compo- 
nents was ejected (jTavlor et al.lll986T) . This event al- 
lowed an accurate measurement of t he jet velocity nea r 
1500 km s" 1 . In HST observations ([Eyres et all (20021 ) . 
arcs can be detected that also could be produced by 
episodic ejection eve nts. X-ray emission was first de- 
tected by EXOSAT (jLeahv fe Tavlorill987l ). and subse- 
quent ASCA observations revealed a complex X-ray spec- 
trum with two soft components (kT = 0.2 and 0.7 keV) 
associated with the jet, an absorbed hard co mponent (7.3 
keV) and a Fe Ka line (jEzuka et all 1 19981 ). They inter- 
preted the hard component as thermal emission by ma- 
terial being accreted on to the white dwarf and the soft 
component as either coronal emission from the giant star 
or emissions from shock s in the jets. Analysis o f arch ival 
Chandra ACIS data bv lGallowav fc Sokoloskil (pool re- 
vealed faint extended emission to the south, aligned with 
the optical and radio jets seen in HST and VLA ob- 
servations. IWheatlev fc Kallmanl (|2006f) reanalyzed the 
ASCA data and interpreted the soft emission as scatter- 
ing of the hard component in a photo-ionized medium 
surrounding the white dwarf. They claim that no other 
sources than the accreting white dwarf are required to 
explain the spectrum. The obvious existence of the jet, 
however, and furthermore the apparent decline of the 
hard X-ray component ob served with the US -Japanese 
X-ray satellite Suzaku by iMukai et al.l (|2006f ) together 
with the lack of a corresponding decline in the soft com- 
ponent, suggest that this interpretation is implausible. 
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Recently, iKarovska et al.l (|2007l . hereafter KCR07) re- 
ported the detection of multiple components, including 
an arc, in the archival Chandra images. 

R Aqr and CH Cyg are the only two jets of symbiotic 
stars which are detected in X-rays. While these two ob- 
jects are seen at high inclinations, in MWC 560 the jet 
axis is practically parallel to the line of sight. This spe- 
cial orientation allows one to observe the outflowing gas 
as line absorption in the source spectrum. With such 
observations the radial velocity and the column density 
of the outflowing jet gas close to the source has been 
investigated in great detail. In particular the acceler- 
ation and evolution of individual outflow components 
and jet pulses has been probed wi th spectroscop i c mon - 
itoring programs, as described in ISchmid et al.l (|2001f) . 
Using this optical data, sophisticated numerical mod- 
els of this pulsed propagating jet have be e n dev eloped 
(|Stute. Camenzind fe Schrniof 12001 IStutel [200l here- 
after Paper I and II in this series) . A number of hydro- 
dynamical simulations (with and without cooling) were 
made in which the jet density and velocity during the 
pulses were varied. The basic model absorption line 
profiles in MWC 560 as well as the mean velocity and 
velocity-width are in good agreement with the observa- 
tions. The evolution of the time-varying high velocity 
absorption line-components is also well modeled. These 
models not only fit the MWC 560 data, but are also able 
to explain properties of jets in other symbiotic systems 
such as the observed velocity and temperature of the CH 
Cyg jet. 

So far, MWC 560 has not been detected in X-rays 
(jMurset et al.lll997l ). We find using the PIMMS tool 1 
that the non-detection in the ROSAT all-sky survey sets 
an u pper limit of the ab sorbed X-ray flux of 0.07 counts 
s _1 (|Murset et ail 1 19971 ) and 7 x 10~ 13 erg s" 1 cm~ 2 , 
respectively. 

The jets in all three symbiotic stars show evidence 
of episodicity. Such episodicity in the ejection process 
has been seen in numerical models of the interaction of 
the stellar magneto s phere and the a ccretion disk (e.g. 
iGoodson et al.lll997HMatt et al.ll2002D. In the disk- wind 
scenario fe.g. iBlandford fc Pavnd 119821 : lAnderson et~aTl 
2003) the time-dependent emission is created by time 
variations in the accretion rate of the underlying disk. 
Unfortunately, so far no hydrodynamical models exist 
for explaining the X-ray emission from symbiotic stars. 
As a first step we have therefore used our existing sim- 
ulations, which fit MWC 560, for understanding the ob- 
served X-ray emission properties of MWC 560, CH Cyg 
and R Aqr. 

In Sj2j we briefly describe the numerical models we have 
used. The total X-ray luminosity and its time depen- 
dence is examined in Sj3j After the resulting spectra are 
calculated in SJU we show emission maps in <j5] and ap- 
ply our main results to X-ray observations of CH Cyg, 
MWC 560 and R Aqr in |j6j Finally our conclusions are 
given in SJ7] 

2. THE NUMERICAL MODELS 
2.1. The hydrodynamical simulations 

1 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 



We solved the equations of ideal hydrodynamics with 
an additional cooling term in the energy equation 

§f + V(pv)=0 

d ( p v) 

tt +V(pv®v) = -Vp-pV$ 
at 

+ V(ev) = -pVv-n 2 A(T) 

at 

p=( 7 -l)e. (1) 

with the code NIRVA NA_CP which w as writ t en by 
iZieeler fc Yorkel (|l99l and modified by iThield (p000h 
to calculate radiative losses due to non-equilibrium cool- 
ing by line emission, p is the gas density, p the pressure, 
e the internal energy density, $ the gravitational poten- 
tial, v the velocity and 7 the ratio of the specific heats 
at constant pressure and volume which is set to 7 — 5/3. 
The general capabilities of the code have been described 
in detail in Paper I, for our approximations and assump- 
tions related to the cooling treatment we refer the reader 
to Paper II. We briefly describe the geometry which we 
have adopted in our simulations below. 

We use a cylindrical coordinate system where the jet 
axis corresponds to the z axis and both components of 
the binary system are located in the plane perpendicular 
to this axis. The hot component is located at the origin of 
the coordinate frame; with a binary separation of 4 AU, 
a red giant is implemented. The red giant is surrounded 
by a stellar wind with constant velocity of 10 km s _1 and 
a mass loss rate of 10~ 6 M Q yr _1 . 

The jet is produced within a thin jet nozzle with a 
radius of 1 AU. The initial velocity of the steady jet is 
chosen to 1000 km s _1 and its density is set to 8.4 x 10~ 18 
g cm~ 3 (equal to a hydrogen number density of 5 x 10 6 
cm~ 3 ). These parameters lead to i) a mass loss rate 
of w 10~ 8 M Q yr _1 in the steady jet, and ii) a density 
contrast between the steady jet and the ambient medium 
77 of 5 x 10 -3 and a Mach number of « 60 in the jet nozzle 
at the origin of the coordinate system. Repeatedly each 
seventh day, the velocity and density values in the nozzle 
are changed to simulate the jet pulses which are seen 
in the observations of MWC 560. The duration of each 
pulse is one day. 

Two models (iv' and i') out of our existing set of eight 
models were chosen for computing X-ray emission prop- 
erties. These models represent maximum (model iv') 
and minimum (model i') values of the jet density in the 
pulses. In model iv' (model i'), the jet density in the 
pulses is higher (lower) than the jet density in the steady 
jet. Although model iv' provided the best fit for the op- 
tical data for MWC 560, our work in this paper shows 
that model i' results in X-ray properties which are more 
appropriate for CH Cyg. For both models we used an 
approximate treatment of radiative cooling. The model 
properties including the velocities and densities of the jet 
pulses are given in Table [TJ 

Maps of logarithm of density and temperature for 
model i' on day 162 and for model iv' on day 115 are 
given in Fig. [TJ Both are the last time-steps calculated. 

2.2. Calculating the X-ray properties 

Wc determined the expected X-ray flux using the den- 
sity and temperature maps from the hydrodynamical 
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TABLE 1 

Parameters of the jet pulses 



model 


"pulse [cm 3 ] 


"pulse [km s x ] 


Mjp [Mq yr" 1 ] 


L jet [erg s 1 ] 


i' 


1.25 x 10 6 


2000 


4.66 x 10" 9 


5.88 x 10 33 


iv' 


1 x 10 7 


2000 


3.73 x 10" 8 


4.70 x 10 34 



Note. — The values of the steady jet emission are n = 5xl0 6 cm -3 , v = 1000 kms" 1 , Mj s = 9.33 xlO~ 9 M yr -1 and L jct = 2.93 xlO 33 
erg s — 1 , respectively. Mj s and Mj p are the mass outflow rates of the jet in the steady state and during the pulse, respectively. The duration 
of each pulse is one day, their period is seven days. 




-20.87 -16.93 -13.00 -20.32 -16.66 -13.00 




1.68 4.40 7.13 1.68 4.36 7.04 




z ( AU ) z ( AU ) 



Fig. 1. — Logarithm of density (top) and temperature (bottom) for model i' on day 162 (left) and for model iv' on day 115 (right), shown 
as contour plots and slices along the jet axis (solid) and parallel to the jet axis at r = 2 AU (dashed). 
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Fig. 2. — X-ray luminosity as a function of the evolution time 
of the jet in model i' (solid) and model iv' (dashed) in the energy 
range 0.15 — 15 keV. Most of the energy is emitted below 2 keV, 
i.e. the plots showing the luminosity in the energy range 0.15-2 
keV are similar. 

simulations as follows. We used the atomic database 
ATOMDB with IDL including the Astrophysical Plasma 
Emission Database (APED) and the spectral models 
output f rom the Astrophy sical Plasma Emission Code 
(APEC. [Smith et al.1l2001l) to calculate the emissivity. 
The default abundances in ATOMDB, i.e. 14 elements 
(H, He, C, N, O, Ne, Mg, Al, Si, S, Ar, Ca, Fe , Ni) 
with solar abundances of lAnders fe Grevessel (|1989t ), are 
used. The energy range is divided into bins of 0.01 keV. 
We compute the spectrum and the total flux in X-rays as 
a function of evolutionary time for each of our models. 

We calculate the X-ray emission in the range between 
0.15 - 15 keV, which is exactly the energy range covered 
by EPIC on XMM-Newton and includes that of the ACIS 
instrument (0.2 - 10 keV) and of HETG (0.4 - 10 keV) 
on Chandra. The emission from gas with a temperature 
lower than < 10 6 K is only marginal in this energy range. 

3. THE TOTAL X-RAY LUMINOSITY AND ITS TIME 
DEPENDENCE 

As expected, the high temperatures, created by the in- 
teraction of the jet pulses with previously ejected matter, 
lead to substantial X-ray emission (Fig. [2]). The X-ray 
luminosity in model iv' is higher than in model i'. This 
is a result of the higher density in the pulses and thus 
of the higher kinetic luminosity L jct = 1/2 M v 2 pumped 
into the jet. Furthermore, in model i' about 5% of the 
average kinetic luminosity is radiated in X-rays, but in 
model iv' about 19%. Since the X-rays are emitted by 
shocked material from the fast moving pulses and since 
the X-ray luminosity is proportional to p 2 , compared to 
Lj et being proportional to p, the ratio of the X-ray lu- 
minosity to the kinetic luminosity is proportional to p. 
Therefore this ratio is higher in model iv' than in model 
i'. 

We find minima and maxima in the X-ray emission Lx 
(computed by integrating over the energy range 0.15 - 15 
keV) which are connected with the periodic emergence 
of jet pulses (Fig. [3]). Thus the period of the variations 
in the X-ray emission is about 7 days. The size of the 
fluctuations is 50 % and more of the average emission. 

While the X-ray luminosity stays constant with time 
for model i', it decreases with time for model iv'. This 
difference might be related to a larger amount of cooling 
in model iv'. The initial shock temperature is identical 
in both models, since the velocities are the same. The 
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Fig. 3. — X-ray luminosity as a function of time for model i' 
(top) and model iv' (bottom); one can see minima and maxima 
created by the emergence of each new jet pulse; the dotted lines 
in the top plots show the different trends of the X-ray luminosity 
decreases with time. 

higher densities in the jet pulses in model iv', however, 
lead to higher densities in the X-ray emitting material 
and thus to higher pressures which result in stronger adi- 
abatic expansion and hence enhanced adiabatic cooling. 
Radiative cooling is also enhanced by the higher densities 
in model iv'. 

4. THE SPECTRUM AND ITS TIME DEPENDENCE 

The spectra of both models in the energy range be- 
tween 0.15 - 15 keV show many different features. They 
show continuum emission, and superimposed on the con- 
tinuum, a large number of emission features (some of 
which are blends of several emission lines). A prominent 
feature, which is mainly due to blended iron lines, is seen 
between 0.7-1 keV. Iron also produces a strong emission 
feature in the 6.4 - 6.7 keV range requiring very high 
temperatures (~ 10 7 K) that are reached locally in the 
jet. 

Like the total X-ray luminosity, the spectrum is also 
highly time-dependent (Fig. 0J. We define two proxies 
for the temperature, one using the low energy spectrum 
and one using the high energy spectrum. These proxies 
can then be used conveniently for direct comparison with 
the single-temperature thermal plasma models typically 
used to fit the observed data. 

4.1. Defining temperature proxies 

In order to characterize the temperature of the prop- 
agating jet from the low energy spectrum, we use the 
fact that below energies of about 0.7 keV, both spectra 
in Fig. 2] are almost identical, but differ significantly be- 
tween 0.7 and 2 keV. Therefore we define the proxy £ for 
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Fig. 4. — Spectrum in the energy range between 0.15 - 15 keV; 
top: for model i' on days 153 (a minimum of the total X-ray lu- 
minosity, solid) and 155 (close to the maximum X-ray luminos- 
ity, dashed), bottom: for model iv' on days 105 (solid) and 107 
(dashed). The spectrum plotted with a solid line is shifted down- 
wards by a factor of 100 for clarity in each plot. One can clearly 
see the time dependence of the spectrum. 

relatively low temperature plasma (10 7 K) in the jet as 



c = 



0.2-0.7) keV 



/(0. 



7-2) keV 



(2) 



with /( .2-o.7)kcV and /( .7-2)keV the fluxes integrated 
over the given energy ranges (Fig. O top). 

The slope of the continuum also changes with time. To 
measure this behavior, we define a flux ratio 



i = 



4kcV 



9kcV 



(3) 



We choose regions in the spectrum where no lines are 
present, although photons with these high energies have 
not been observed from the jet or might be confused with 
photons from the central engine in observed spectra. 

4.2. Determining the temperature from model spectra 

In model i', the high temperature proxy, log(£), varies 
with time spanning a range from 1.5 to 3. It has its high- 
est value (i.e. the spectrum has the steepest slope, hence 
the average temperature is at its lowest value), when the 
X-ray luminosity shows a minimum (Fig. [6J. Compar- 
ison of the model log(£) values with that of a single- 
temperature thermal plasma (Fig. El bottom) gives us 
temperature estimates of the hot component between 
8 x 10 6 K (0.69 keV) and 1.7 x 10 7 K (1.5 keV). The 
highest temperatures are consistent with that of post- 
shock gas with a shock velocity of about 1100 km s _1 . 
The minimum in the log(£) (i.e. maximum in average 




T ( K ) 

Fig. 5. — Temperature proxies as a function of temperature for 
a single-temperature thermal plasma calculated with ATOMDB: 
hardness ratio £ (top) and flux ratio £ (bottom) as defined in the 
text. 




Fig. 6. — X-ray luminosity of the jet as a function of time in the 
energy range 0.15 - 15 keV (top), flux ratio £ as a function of time 
(middle) and hardness ratio £ as a function of time (bottom) for 
models i'. 

temperature) coincides with the emergence of each new 
pulse; within the next 2-3 day period the compressed 
knot cools and the emissivity increases. Therefore the 
maximum in the X-ray luminosity is reached about 2-3 
days later. 

The low temperature proxy, log(£), varies between 
about 0.5 and 1.6, the corresponding temperatures lie 
between 1.6 x 10 6 K (0.14 keV) and 3 x 10 6 K (0.26 
keV). The jet is therefore better described as a combi- 
nation of a warm and a hot component rather than as a 
single-temperature plasma. 

In model iv', log(£) varies between 2.1 and 2.8; the 
corresponding temperatures are 8 x 10 6 K (0.69 keV) 
and 1.2 x 10 7 K (1.03 keV). The low temperature proxy, 
log(C), lies between 0.1 and 0.6; the corresponding tem- 
peratures are 3 x 10 6 K (0.26 keV) and 3.8 x 10 6 K (0.33 
keV). As in model i', the jet is better characterized as a 
combination of a warm and a hot component. 

The range of temperatures in the hot component over 
one pulse cycle in model iv' is smaller compared to that 
in model i'; this is because the higher density in model 
iv' makes radiative cooling more efficient, such that the 
shock heating is damped more efficiently. The different 
density contrasts between the jet pulses and the steady 
jet in both models also lead to different shock velocities 
and thus to different shock temperatures to which the 
plasma is heated initially. 
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Fig. 7. — same as Fig. \E\ but for model iv' 
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Fig. 8. — Spectrum in the energy range between 0.15 - 15 keV for 
model i' on days 155; several strong emission lines are present with 
very different ionization potentials which can only be explained 
using two components at different temperatures. 

4.3. Emission lines in the model spectra 

We have identified the strongest emission lines char- 
acterizing a typical jet spectrum using the spectrum de- 
rived from model i' on day 155 as a template (Fig. [8|). 
The most prominent lines are the oxygen line at 0.57 
keV, the Ne lines at 0.93 and 1.03 keV, the Mg lines at 
1.35 and 1.47 keV, the Si lines at 1.86, 2.01 and 2.18 
keV, and the Fe complex at about 6.5 keV. All strong 
lines with their identifications are given in Table in- 
cluding their fluxes at days 153 and 155 in model i' and 
days 105 and 107 in model iv'. Most of these lines are 
hydrogenic or He-like lines of heavy elements, however, 
also lines of higher ionization states are present (Fe XXII, 
Fe XXIII, Fe XXV). Thus this set of lines with very dif- 
ferent ionization potentials can also only explained with 
two components at different temperatures. 

5. THE X-RAY EMISSION MAPS 

The results in the previous section suggest the exis- 
tence of two components, a warm one with temperatures 
in the range of (1.6 - 3.8) xlO 6 K and a hot one with 
temperatures of (8 - 17) x 10 6 K, respectively. As already 
pointed out in Paper I and II, the jet consists of dense, 
cool knots and tenuous, hot inter-knot gas. The knots 
are too cold to emit X-rays, thus the low and high energy 
components in the X-ray spectrum probe the tempera- 
ture structure of the hot parts of the jet (Fig.l and Fig. 
12 in paper I). These hot parts consist of shocked mate- 
rial in the inter-knot gas segments. Temperature gradi- 



ents are present within each inter-knot gas segment and 
also between the inter- knot gas close to the jet source and 
those downstream near the jet head. The regions with 
the highest temperatures and thus emission above 6 keV 
lie in the first two inter-knot gas segments within only 
a few AU from the central source (Fig. [9]). The inter- 
knot gas segments become progressively cooler, as they 
move farther away from the jet source due to adiabatic 
expansion and the resulting cooling of jet material. 

6. COMPARISON WITH OBSERVATIONS 

We now compare the results of our models to X-ray 
observations of the three objects MWC 560 (only upper 
limits to the flux are available), CH Cyg and R Aqr. 

6.1. MWC560 

The source fluxes in the 0.2 - 2.4 keV range are 
3 x 10 -13 erg s" 1 cm~ 2 for model i' and 2 x 10~ 12 
erg s _1 cm~ 2 for model iv', us i ng a distance of 2.5 
kpc to MWC 560 (jSchmid et all 120011 . and references 
therein). The latter flux is reduced to an absorbed flux 
of 1.7 x 10 -13 erg s" 1 cm~ 2 , using the mod el of the visual 
interst ellar extinction in the Galaxy by Hakkila et al. 
(fl997l ) which gives an A v = 0.88 or N H = 1.55 x 10 21 
cm~ 2 for MWC 560. The absorbed flux is consistent with 
MWC560's non-detection in the ROSAT all-sky survey 
(< 7 x 10~ 13 erg s" 1 citT 2 , iMiirset et a!]|1997D . The 
fluxes, however, are high enough such that today's X-ray 
telescope as Chandra and XMM-Newton should be able 
to detect MWC560 2 . 

6.2. CH Cyg 

Since the jet velocities are similar in MWC 560 and CH 
Cyg, we can make a detailed comparison of our models 
with the latter object. This comparison is limited to the 
soft emission from the propagating jet and excludes the 
X-ray emission above about 2 keV which is believed to 
be dominated by the variable scattered hard X-rays from 
t he central source. 

Ezu ka et al.l (|1998h resolved for the first time atomic 
emission lines (or their blends) from elements in hydro- 
genic and He-like ionization states in the X-ray spectra 
of CH Cyg. The most prominent features in the observed 
spectra are the oxygen line at 0.57 keV, the blend of i) 
Ne lines at 0.93 and 1.03 keV, ii) Mg lines at 1.35 and 
1.47 keV, iii) Si lines at 1.86, 2.01 and 2.18 keV, and the 
Fe line complex at about 6.5 keV. These lines are also 
seen in our model spectra. Furthermore, they also had 
to introduce a two-temperature thermal plasma model 
to explain the set of lines detected in the ASCA spec- 
trum. The temperatures of their two components (0.21 
and 0.72 keV) are in the same range as in our model i' 
(warm component: 0.14 - 0.26 keV, hot component: 0.69 
-1.5 keV, see §4. 2ft . In model iv', the warm component 
is too hot (0.26 - 0.33 keV) to explain the observations. 

Recently, KCR07 reported the detection of multiple 
components, including an arc, in the archival Chandra 
images. This arc has a similar opening angle to many 
of the arcs visible in our X-ray emission maps (Fig. [9]), 
and its presence supports our pulsed-jet model in which 

2 XMM-Newton observations proposed by the authors will be 
executed in AO-6 after May 2007. 
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TABLE 2 

Prominent emission lines in the model spectra 



ion energy (kcV) flux (erg s 1 cm 2 ) 

model i' model iv' 
day 153 day 155 day 105 day 107 



C V 


0.31 


8.81x10" 


-14 


7.02x10" 


-14 


1.36x10" 


-13 


7.03x10" 


-14 


C VI 


0.37 


3.09x10" 


-14 


2.49x10" 


-14 


1.06x10" 


-13 


4.21x10" 


-14 


N VI 


0.42 


6.96x10" 


-15 


6.53x10" 


-15 


2.58x10" 


-14 


1.08x10" 


-14 


N VI 


0.43 


1.00x10" 


-14 


9.20x10" 


-15 


4.06x10" 


-14 


1.50x10" 


-14 
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Fig. 9. — Emission maps of model i' on day 155 in the 0.15 - 1.5 keV range, in the 1.5 — 3 keV range, in the 3-6 keV range and in the 
6-7 keV range. 



such emission results from internal shocks generated by 
colliding jet ejecta. 

Considering the smaller distance of 268 pc to CH Cyg, 
the model X-ray luminosity of 2.2 x 10 32 erg s' 1 for model 
i' and 1.5 x 10 33 erg s _1 for model iv' in the 0.2 - 2 keV 
range implies fluxes of 2.6 x 10~ n erg s _1 cm™ 2 and 
1.74 x 10 -10 erg s _1 cm -2 , respectively. The interstellar 
extinction for CH Cyg is A v = 0.006 or JV H = 1 x 10 20 
cm~ 2 , thus the absorbed flux is 2 x 10 -11 erg s _1 cm -2 
for model i' and 1.3 x 10~ 10 erg s -1 cm -2 for model 
iv', respectively. The measured fluxes of the soft compo- 
nents associated with the jet, howev er, lie in the range 
(0.4 - 3.8) x 10~ 12 erg s' 1 cm~ 2 (lEzuka et alJ fl998l 
iGallowav k Sokoloskil [200l iMukai et al.l 12006ft . Since 
these measured fluxes are smaller than those predicted 
by model i' and by far smaller than those from model iv' 
which has a higher jet pulse density than model i', we 
infer that the jet pulse densities in CH Cyg are smaller 
than those in MWC 560 and our models. If, as stated in 
121 the ratio of the X-ray luminosity to the kinetic en- 
ergy pumped into the jet is proportional to the jet pulse 
density, we can estimate that the jet pulse density has to 
be reduced to about 10 5 cm~ 3 to model the X-ray flux 
observed in CH Cyg. In order to resolve the difference 
between measured and model fluxes with uncertainties 
in the distance to CH Cyg, the latter would have to be 
about 700 pc, however, this large value is very unlikely, 
since the distance was measured with HIPPARCOS with 
an error of 23 %. 

KCR07 give an observational estimate of the density 
in the jet, 50 cm~ 3 , based on the total X-ray luminosity, 
assuming an emitting sphere with a radius of 400 AU 
and a mean emissivity of 2 x 10 -23 erg cm 3 s -1 . How- 
ever, astrophysical jets, typically, are collimated struc- 
tures which have significantly smaller volumes than a 
sphere. Since we simulated the propagation of the model 
jets only up to a length of 65 AU, not 400 AU, we scale 
up its volume as follows. If we assume that the cross- 
sectional area stays constant as the jet propagates and 
evolves with time, the jet volume is about (tt x 6 2 x 400) 
AU 3 (Fig. [9]) which is smaller by about a factor of 6000 
than that of a sphere with a radius of 400 AU. Alter- 
nately if we assume that the lateral expansion of the jet 
is proportional to its axial expansion 3 , the volume is still 
overestimated by a factor of about 160. Furthermore, 
within the volume of the jet, the X-rays are emitted by 
clumps and not by the whole jet, i.e. the filling factor 
for the X-ray emission is, /i < 1. In addition, KCR07 
assume a temperature of 1.2 x 10 7 K (1 keV) to estimate 
the emissivity, however, our jet models show that such 
high temperatures are only archieved in the innermost 
region of the jet; at larger distances from the jet source 
along the axis, the X-ray emitting knots are significantly 
cooler (Fig. Q]). Thus the value of the emissitivity should 
be lower and the density should be higher than estimated 
by KCR07 by another correction factor fa > 1. Hence 
we conclude that an accurate estimate of the X-ray emit- 
ting volume and the temperature would lead to a higher 
density by a factor of (160-6000)/2//i compared to the 
value given by KCR07, of the order of 10 4 -10 6 cm -3 . 
This range compares well with the densities of material 
emitting soft X-rays in our models which are of the order 

3 which is seen in our models after day 70 (Paper II) 



of 10 5 -10 6 cm" 3 (Fig. [J). 

Other possibilities to reduce the X-ray fluxes in our 
models, bringing them closer to the observed ones in CH 
Cyg, are i) a longer timescale between the pulses and ii) a 
smaller velocity difference between the steady jet and the 
jet pulses. In the first case, less energy is pumped into 
the jet and each jet knot and the jet head can cool further 
before being hit by the next pulse. In the second case, 
the smaller velocity difference reduces the temperature 
to which the shocked material is heated. However, the 
density contrast then also has to be adjusted in order to 
reproduce the observed proper motion of the jet knots. 
Which of the above scenarios is the most likely one can 
only be determined from future simulations which have 
been fine-tuned to fit the properties of the X-ray emitting 
material in the jet of CH Cyg. 

In 1994 and 2006, the flux from the jet in CH Cyg was 
almost at the sa me level, in the ran g e (2.7 — 3.8) x 10~ 1 2 
erg s" 1 cm" 2 (jEzuka et all 119981 : IMukai et all l2006h . 
In 2001, however, it was lower by a factor of ten 
(|Gallowav k So koloski 2004) . In the context of our mod- 
els, such a drop in flux may result from a large decrease 
in the density of the pulses which may be caused by a 
drop in the accretion rate onto the white dwarf. 

6.3. R Aqr 

Since the jet velocity in R Aqr is smaller by a factor of 
2 than those in our models, we can only make a more 
limited co mparison of the mod el results with the ob- 
servat ions. iKellogg et al.l (|2006f) and also iKorreck et al.l 
(2006) report a tange ntial motion of an X -ray emitting 
knot 4 of 600 km s _1 . IKorreck et al.l (|2006ft estimated a 
density of 100 cm' 3 , however, it may be possible that, 
as in the case of CH Cyg, the X-ray emitting volume is 
o verestimated and thu s the density is underestimated. 

IKellogg et ail (|200lD found in R Aqr that the NE jet 
is more luminous by a factor of 3 than the SW jet. The 
spectrum the NE jet was fitted with a single-temperature 
thermal plasma with a temperature of 1.66 keV, the spec- 
trum of the SW jet with a plasma temperature of 0.2 
keV. The simplest explanation of this difference is that 
the ambient media on both sides have different densities 
which would lead to different compression factors and 
thus to different shock heating temperatures. However, 
proper motion measurements show no significant differ- 
enc es in the velocities of the knots in both side s of the 
jet (jParesce k Hackl 119941: iMakinen et al.ll2004D . which 
we would expect if different ambient densities would de- 
celerate the jet differently on both sides. 

Our modeling can provide a plausible explanation for 
the observed differences between the two sides of the jet, 
if we assume that the ejection of the jet pulses on both 
sides are out of phase with each other. iHollis et al.l (|199lD 
derived a kinematic age of both jets of about 90 yrs; over 
the extent of the jet we can observe 3-6 knots in the ob- 
servations of e.g. iParesce k Hack! (|1994ft which suggests 
several ejection events in this period of time. A period 
of about 17 years for these ej ection events has be en in- 
ferred from radio observations (|Kafatos et al.lll989l ). thus 
they are significantly larger than in MWC 560 and in our 

4 IKorreck etHTl ^OOS ) give an even smaller velocity of 380 km 
s _1 in their Table 1, however, without explaining the difference 
from the value given in the text. 
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models 5 . We hypothesize that the X-ray emitting blobs 
in the NE jet were ejected later than those in the SW jet 
and therefore have cooled less. This hypothesis is sup- 
ported by the fact that the X-ray emitting component in 
the NE jet is closer to the central core (about 8" ) tha n 
the blobs in the SW jet ( 1 2-26" . [Kellogg et al.ll200l . 
A new SW jet component with an offset of about 1.5" 
from the cen t ral so urce has recently been reported by 
iNichols et al.1 (l2007f). Assuming a jet velocity of about 
600 km s" 1 (|Kellogg et al.l I2006D . i.e. about 0.6" per 
year, we obtain a kinematic age of about 2.5 years for 
this component. Even if t he new component h ad been 
ejected during the epoch of lKellogg et al.1 (|2001h 's obser- 
vations, its flux would have contributed to the central 
source, but not to that of the SW jet. Hence the pres- 
ence of the new component does not conflict with our 
hypothesis. The time period between the emergence of 
the new SW jet component and the emergence of the blob 
now located at 12" is about 17 years which supports the 
inferred period for the ejection of jet pulses in R Aqr. 

In we showed that the total X-ray luminosity de- 
creases with time, probably due to adiabatic cooling in 
the jet (Paper II). This effect provides a plausible ex- 
planation of the decrease in X-ray flux in the jet of R 
Aqr fr om a value of 5 x 10 ~ 13 erg s _1 cm -2 in the early 
1990s (Hi insch et aT1ll998h to 1 x 10" 13 erg s" 1 cm" 2 in 
2000 (jKellogg et al.ll2001h . If this explanation holds, we 
expect to see a further decrease in flux in future obser- 
vations. 

6.4. The 6.4 ~ 6.7 keV iron line complex 

This iron line complex has been observed in both ob- 
jects, CH Cyg and R Aqr. Our model spect ra also show 
the ex istence of this Fe line complex (Q. lEzuka et al.l 
(1998) fitted the observed spectrum of CH Cyg with three 
single-temperature components (a warm and a hot com- 
ponent to explain the jet emission and hard component 
for the central engine) and an additional Gaussian rep- 
resenting fluorescence emission in the Fe Ka line. This 
fluorescence occurs close to the white dwarf and the ac- 
cretion disk. Since it is not possible to disentangle the 
flux of the thermal and the fluorescence components in 
this line complex and since our models do not include 
the effect of fluorescence, we cannot compare our models 
with this part of the observed spectrum of CH Cyg. 

In R Aqr, the origin of the hard X-ray emission is more 
ambiguous. ;Kello gg et al.l (|2001h detected 16 photons in 
the range between 6.4 and 6.7 keV which they attribute 
to the central source due to the extraction regions they 
chose. The physical origin of this emission, i.e. ther- 
mal or fluorescence, cannot be decided, since there are 
not enough photons to model the spectrum in this en- 
ergy and characterize its nature. In our model, we find 
significant hard radiation including continuum and iron 
emission lines, being emitted by the first two internal 
shocks in the jet downstream from the source (i.e. at a 
distance less than 15 AU). Since at R Aqr's distance of 
200 pc, 15 AU correspond to 75 milliarcseconds, which is 
well below the angular resolution of Chandra, the X-ray 
emission from these shocks cannot be separated from the 

5 The events in R Aqr are thought to be triggered by periastron 
passage, while the variations in MWC 560 seem to be a result of 
disk instabilities. 



central source. The model X-ray luminosity in the 6.4- 
6.7 keV range is between 5.1 x 10 24 erg s _1 and 7.3 x 10 28 
erg s _1 for model i' and between 2.8 x 10 28 erg s _1 and 
5.1 x 10 29 erg s _1 for model iv', respectively, depend- 
ing on whether the jet is in its minimum or maximum 
state. At a distance of 200 pc, this corresponds to fluxes 
between 1.1 x 1 0~ 18 erg s -1 cm" 2 a nd 1.1 x 10~ 13 erg 
s _1 cm -2 . Since iKellogg et a l. (2001) measured a flux of 
4.9 x 10~ 14 erg s _1 cm -2 at 6.41 keV, we suggest that 
the measured iron line flux may be emitted by the jet 
itself, and an additional fluorescence component is not 
needed, in R Aqr. However, new simulations with the 
same jet velocity as in R Aqr are needed in order to test 
this suggestion. 

7. CONCLUSIONS 

We have used our models of pulsed, radiative jets in 
symbiotic stars in order to investigate their X-ray prop- 
erties in detail. These models show that the well-studied 
pole-on jet in MWC 560 should be easily detected by 
today's X-ray telescopes such as Chandra and XMM- 
Newton, since our model flux and its time variation for 
this source are of the order of 10 -13 erg s _1 cm~ 2 . 

We find minima and maxima in the X-ray emission Lx 
(computed by integrating over the energy range 0.15 - 15 
keV) which are connected with the periodic emergence 
of jet pulses. The maxima of the total X-ray luminos- 
ity occur 2-3 days after the emergence of new jet pulses, 
which are ejected every 7 days. The size of the fluctua- 
tions is 50 % and more of the average emission, making 
such X-ray flashing jets detectable with Chandra and 
XMM-Newton. 

The X-ray spectra of our model jets are rich in emission 
line features, the most prominent of which correspond to 
observed features in the spectra of CH Cyg. 

By using low and high energy temperature proxies de- 
rived from the spectra, we can show that the emission 
can be adequately characterized with a hot and a warm 
optically-thin plasma component. The hot component 
has temperature values of about 0.7 keV (1.6 keV) during 
the minima (maxima) of the total X-ray luminosity and 
the warm component has temperature values of about 
0.14 keV (0.33 keV) during the minima (maxima). 

While model iv' is appropriate for MWC 560, we have 
shown that model i', which has a lower jet pulse density 
than model iv', is more appropriate for the jet in CH Cyg 
in terms of explaining the lower X-ray flux. Other possi- 
bilities to reduce the flux are a longer timescale between 
the pulses and a smaller velocity difference between the 
steady jet and the jet pulses. Which of the above sce- 
narios is the most likely one has to be tested in future 
simulations which have to be tuned to the jet in CH Cyg. 

Our models also offer a plausible explanation for the 
differences in luminosities and temperatures in the NE 
and the SW jet of R Aqr. We assume that the ejection 
of the jet pulses on both sides are out of phase with each 
other. We hypothesize that the X-ray emitting blobs in 
the NE jet were ejected later than those in the SW jet 
and therefore have cooled less. 

We find the existence of iron line emission in the 6.4 
-6.7 keV range in our models which is also observed in 
both, CH Cyg and R Aqr. Our models cannot be directly 
applied to CH Cyg, because of an additional fluorescence 
component from the central source and accretion disk in 
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the observed spectrum and because fluorescence is not in- 
cluded in our models. In the case of R Aqr, although this 
emission has been associated with the central source, our 
modeling shows that it is consistent with being produced 
by jet knots very close to the latter, because their sepa- 
ration in our model is well below the angular resolution 
of Chandra. 

Using our current models which were built to fit the 
optical data of the jet in MWC 560 we are able to explain 
some of the important characteristics of X-ray emission 
from jets in MWC 560 and other symbiotic stars. The 
results of this study demonstrate the great potential of 
future numerical simulations of pulsed jets which have 
been fine-tuned to specific source properties for under- 
standing the jet phenomenon in symbiotic stars. Fur- 



thermore, we hope that this study will inspire new and 
more sensitive observations of X-ray emission from jets 
in symbiotic stars. 
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